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Abstract. 

We study the possibility of producing and detecting continuous variable cluster 
states in an optical set-up in an extremely compact fashion. This method is based 
on a multi-pixel homodyne detection system recently demonstrated experimentally, 
which includes classical data post-processing. It allows to incorporate the linear optics 
network, usually employed in standard experiments for the production of cluster states, 
in the stage of the measurement. After giving an example of cluster state generation 
by this method, we further study how this procedure can be generalized to perform 
gaussian quantum computation. 
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1. Introduction 

Quantum computing is a very promising subject despite the formidable challenges it 
must overcome [I]. One of those challenges lies in the difficulty of finding a system 
allowing for a large number of entangled states to be manipulated and for many quantum 
gates to be implemented successively. In this context, an interesting avenue consists in 
measurement-based quantum computing [21 [3j H] (MBQC) protocols which rely on the 
availability of a large, multipartite entangled state on which a series of measurements 
is performed. For each operation that we wish to implement, a specific sequence of 
measurements on its modes has to be chosen. The measurement outcomes are in general 
used to determine which measurement has to be performed afterwards, and to perform a 
final correction stage on the output state. Optical systems are promising candidates for 
the experimental realization of MBQC j5j|6] and convincing demonstrations have already 
been made, both in the discrete [7] and continuous variable regime [8]. Continuous 
variable (CV) systems are especially promising since they allow to create and to process 
cluster states in a deterministic fashion, contrarily to the discrete case [9]. 

Most often, the creation of multimode entangled states such as cluster states in 
a CV quantum optical setup is implemented using a series of squeezers followed by a 
network of beam-splitters and dephasers, which transform the squeezed input modes 
into entangled output modes pm El HU [12] ■ The configuration of this network varies 
considerably with the state to be generated, and its complexity grows rapidly with the 
number of modes, which renders this method poorly scalable. 

In a recent experiment at Australian National University [13], some of us 
have experimentally demonstrated the possibility of performing exceptionally compact 
operations on optical modes using spatial (transverse) modes of light. In this experiment, 
multiple transverse modes were generated and individually squeezed in optical cavities, 
and then mixed together to produce a multimode squeezed beam. This beam was then 
measured by means of a multi-pixel homodyne detection (MPHD) system, followed by 
digital post-processing of the acquired signals, e.g. multiplication of each signal a by 
suitable gain [H]. This allowed them to measure simultaneously an arbitrary quadrature 
on all the modes in a chosen mode set. The use of this method to emulate cluster states 
statistics was also sketched in this work. 

It is also possible to generate multimode quantum states of light in the frequency 
domain starting from an optical femtosecond frequency comb and using a synchronously 
pumped optical parametric oscillator (SPOPO). In this context, the formation of three- 
mode squeezed non-classical states has already been demonstrated [15J. A multimode 
detection system analogous to the one employed in Ref. [13] can be implemented for the 
detection of these frequency modes, where the spatial pixels are replaced by frequency 
bins. An appealing advantage of this setup is that the source is now scalable, as all the 
modes are produced within only one optical device, and potentially hundreds of modes 
can be entangled. 

This motivates our study of the ensemble of equivalent unitary operations which 
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can be readily implemented by the multi-pixel detection followed by digital signal 
recombination. In particular, we give a procedure to assess whether and how a given 
cluster state can be realized with this method. Then, we apply this idea to perform 
gaussian quantum computation in the measurement-based model and give some specific 
examples. 



2. Modelization of the operations which can be performed on the modes 

Let us first model an experiment that would consist in generating a multimode squeezed 
state and in performing a multi-pixel homodyne detection from the point of view of 
canonical transformations on the optical input modes. Such an experiment can be 
separated in four steps, as can be seen in Figjlj a) generating a composite beam made 
of independent orthogonal squeezed modes, which can be either spatial or frequency 
modes; b) mixing it on a beam splitter with a shaped local oscillator; c) propagating it 
to multi-pixel detectors; d) recombining digitally the acquired signals. 

As discussed above, there are several ways of generating multimode squeezed beams. 
Nevertheless, independently of the choice adopted, the result can be modeled as an 
optical field composed of iV orthogonal and normalized modes of the electro-magnetic 
field Ui(p), where p could either be spatial coordinates (x,y) or frequency components 
bj depending on the realization. We assume that in each mode the state of the field is 
infinitely squeezed, say along the p quadrature, corresponding hence to the p-eigenstate 
\s) p , being p\s) p = s\s) p , with zero eigenvalue, i.e. |0) p = l/y/2~7i J ds\s) q [16]. The 
discussion of the effect of finite squeezing is beyond the scope of this paper. 

These N modes are accompanied by vacuum modes in all the other modes of the 
complete modal basis {ui(p)}. To each mode are associated a pair of bosonic operators 
d Ui and a) u . with a Ui = J dpd(p)u*(p) (where J dp stands for JJ^ gR2 d 2 p if p is a 2D 
variable), and a(p) = J2T=i a u k Uk(p)- These input modes can be grouped in a vector 
a u = (a u , 0j\j) = (a ui , a,^, d M1 , o U2 ...) • 



Once the multimode squeezed beam is generated, it is mixed on the beam splitter 
with a local oscillator. Then, it is propagated to two arrays of detectors composed of 
pixels, each of surface Si, that collect the intensity of the incoming light from the two 
respective beams, and the differences pixel by pixel of the two signals are taken. This 
propagation is governed by diffraction in the spatial case, and by frequency dispersion 
realized with the help of a prism in the frequency case (see Fig{TJ. To describe the effect 
of the detection of the beam it is convenient to introduce another set of modes, that we 
shall call pixel modes; these are defined as 

v f p \ _ J Ki u lo(p) Vp G Si 

I elsewhere. ' 

and correspondingly d v . = J dpd(p)v*(p), where Ki is a normalizing constant such that 
I P m^P v i^ = K< i ffp^s d 2 P I m lo(p)| 2 — 1- Ec[-([l| implies that each pixel mode is a 
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Figure 1. a) Sketch of the experiment performed with frequency modes. The 
(infinitely) p-squeezcd modes Ui(p) generated by synchronously pumped optical 
parametric oscillator (SPOPO) are mixed with a local oscillator on a 50/50 
beamsplitter. Both outputs are then dispersed and each frequency band is sent to a 
different pixel of the two multi-pixel detectors. The differences of the acquired signals 
are taken pixel by pixel. These are later processed by an ordinary computer which 
multiplies every trace by a suitable gain, for instance to obtain a cluster state. The 
additional correction step indicated in the dashed box is only required when performing 
a quantum computation; it reinterprets the measurement result on the mode of interest 
depending on the outcomes of the measurements on the other modes. Note that an 
analog device can similarly be implemented for the detection of spatial input modes, b) 
Functional modelization of the same experiment. Aopo models squeezing quadratures 
of source modes, Ut light propagation from source to multi pixel detector, Alo local 
oscillator phases in the pixel basis and O digital data processing. Refer to main text 
for details. 

"slice" of the local oscillator, i.e. it shares its shape but only on a window with width 
given by the pixel dimension (and it is zero ouside). Note that the pixel modes form a 
complete basis, allowing to reconstruct all the possible modes, only in the limit of infinite 
pixel number. The detection process can be modeled as an effective transformation 
corresponding to light propagation from the squeezed modes set to the pixel modes 
one, and to a homodyning in that basis with an intense local oscillator «lo,iG°)- Here 
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for consistency we assume that p is, from the beginning, the variable that spatially 
describes the field in the multi-pixel detector plane. For instance, in Fig. [T] it does 
indeed correspond to frequency co, but in the spatial case it depends on the actual 
imaging system between squeezing generation and the detectors. This transformation 
is described by a matrix Ut acting on the bosonic operators: a v = Uto u , which brings 
the input modes a Ui onto the pixel modes a Vi , with Unj = K i f e g. dpu^ Q (p)uj(p). More 



details can be found in Appendix A. 2 and in Refs. [T71 US]. Note that in principle the 
detection matrix Ut does not need to be square, since the number of pixels P may differ 
from the number of input modes - the latter may be potentially infinite. 

As an example to fix the ideas, that we will use later, consider the simple case 
in which four input squeezed modes are detected with the help of four pixel modes. 
Let us assume that the intensity distribution of the local oscillator mode is the first 
mode of this basis, i.e. mlo(p) = u i(p)- We take a toy model for the input modes, 
describing the mode n by a square profile in which we introduce n — 1 7r-phase shifts 
(or "flips") at regular intervals, as illustrated in Figj2] (top panels). These modes are 
however qualitatively similar to the modes generated by a SPOPO, which in turns are 
close to Hermite-Gauss modes [T5]. In this simple case, the matrix Ut is unitary and 
takes the simple form 

1 \ 

-1 



/ 1 



Un 



\ 



1 



1 1 

1 -1 



(2) 



j = o 




i = 



j = l 



j = 2 



j = 3 




i= 1 



i = 2 



i = 3 



Figure 2. Input square flip- modes Uj(p) (top panels) and pixel modes i>j(p) = 
^2iUTijUj(p) (bottom panels) for < i,j < 3, as a function of p. As it can be 
seen, the latter allow to reconstruct the profile of the former and conversely. 



Above we have assumed that the input state is p-squeezed on each mode. 
Experimentally it is in principle possible to modify the squeezing quadrature. For 
simplicity, we model that by changing the input mode phase instead of its state, which 
is completely equivalent in an homodyne detection scheme. This is done by introducing 
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a diagonal dephasing matrix Aopo with Aopoij = e^dij- The input modes entering 
the detection system are u'^p) = e _ *Uj(p) and a u > = e l ^a u . The transformation of 
these input modes to the pixel set as defined above is now a v = UTA* -p a u i , as Aqpq 
moves the phase shifted modes to the modes a u and Ut is the transformation from a u 
to a v . 

With the same approach, the ability of shaping the phase of the local oscillator 
by phase-shifting the pixel modes can be modeled by the action of a diagonal matrix 
a v i = Alo<^ with elements Aloij = em 8i,j, leading to the modes v[(p) = e~ tipi Vi(p). 

Finally, after the detection, it is possible to digitally recombine the electronic signals 
coming from each pixel by multiplying them with real gains [H] . This allows us to 
detect the field according to desired modes (e.g., as we shall see later, according to the 
modes which correspond to a cluster state), and amounts to applying a P-dimensional 
orthogonal matrix, i.e. a v " = Oa v i = a out (see also the Appendix). 

In summary, we find that the series of the possible transformations that we can 
perform on the input modes are 

flout = OAloG a u = ?7mphd 4, (3) 

where A LO = diag(e^ 1 , e lLP2 , . . . e lipN ), O is real orthogonal, and G = U T A* OPO is a fixed 
matrix. Note that since AloA lo = X and GG^ = X the product OA^qG = ?7mphd is 
unitary. As implicit in the previous discussion, the matrices O and A L o are easily 
tunable in the experiments. The dephasing matrix Aqpo can also be adjusted to some 
extent by changing the spectral shape of the pump [18] . 



3. Characterization of the feasible transformations 

We now address the characterization of the sub-set of unitary operations, acting on the 
annihilation operator a u , which can be emulated by means of Eq.(|3]). The question: 
"which class of unitary transformations can we emulate in the laboratory by means 
of the multi-pixel homodyne detection plus data processing?" can be recast into the 
following question: "given a unitary matrix U t h, can it be written in the form of Eq. ([3])?" 
Answering this questions requires finding a solution for the parameters of O and A L o 
(given a certain fixed detection matrix G) such that 

Uth = Umpud- (4) 

It is easily seen that a necessary and sufficient condition for Eq.Q is 

U'JU'^ = D, (5) 

with U' th = UfaG* and where D is a diagonal matrix with unit modulus complex 
elements. The necessity is simply proved by assuming that Eq.([4|) holds, and by 
direct calculation of the product U^U^ upon substitution of Eq7^3|, which yields 
U' t ^U' th = G*G T Al T OA h0 GG ] = A^ Q , where we have used A£ = A LO since this 
matrix is diagonal, and where Alo is obviously a diagonal matrix with unit modulus 
complex elements. For the sufficiency, we can show that if Eq.(|5]) holds, we can always 
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find at least a set of "experimental" matrices Alo and O such that their combination of 
the form of £/mphd in Eq.pl) satisfies Eq.Q. This is achieved by taking as a realization 
of the diagonal matrix Alo any of the solutions of the equation 

A LO = D^ = (Ui T U{ h )l (6) 

Then, it is automatically ensured that, for each solution Alo, the matrix 

O = £4Al (7) 

is orthogonal and fulfills Eq.rtiJ. Indeed by direct calculation we obtain £/mphd = 
OA LO G = L^AlqAloG = U t h, which completes our proof. 

As a major result of our work this shows that, given a certain unitary matrix 
f/th, if a solution of Eq.([5]) exists, with the use of Eqs.([6]) and we can determine 
the suitable experimental parameters in terms of the gain coefficients O and the local 
oscillator phase shaping Alo- These allow to obtain with the multi-pixel homodyne 
detection measurement the statistics of the quadrature outcomes as if the matrix U t h 
- usually corresponding to a complex optical network - had been applied to the input 
squeezed modes, and the same quadratures were measured one by one. Let us also note 
that our proposal is versatile, in the sense that different matrices U t h can be produced 
with little or no reconfiguration of the experimental set-up. 

3.1. Cluster states in the multimode beam 

To provide a demonstration of our method, we consider the simple case in which four 
input squeezed modes are detected with the help of four pixel modes, which entails the 
mode transformation given in Eq.Q. Let us consider as a concrete example the case in 
which the desired unitary transformation to be emulated is the one which transforms N 
squeezed modes in a linear cluster state. 




Figure 3. Four-mode linear cluster state. Each circle represent a cluster mode and 
the edges represent entanglement between modes, namely that the two modes have 
been acted upon with a C^-type interaction e 1 ^®^ [6]. 



In the case N = 4 (Fig. [3]), such a transformation is represented by the following 
matrix [TU] 
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In order to find a simple algebraic solution for the problem we consider the dephasing 
diag(l,M, 1). With the use of Eqs.Q, Q, @ and g we find 



matrix Aq PO 



TV T TV 
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-2-i 










2-i 
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-2+i 
V5 
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J 



(9) 



Since the matrix in Eq.([9]) is diagonal, the condition in Eq.([5]) is satisfied, which implies 
that the corresponding cluster state matrix ^ can be implemented experimentally. 
Among the 2 N = 16 possible solutions for Alquh an d the corresponding orthogonal 
matrix 0\ in satisfying Eqs. (Jsl Pfl) for U t h = Un n we can chose for instance the set 
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(10) 



with 7 = arctan | and a = cos 1 ^ — ^gj. Hence, we conclude that shaping the 

phase of the local oscillator on each pixel modes as prescribed by the matrix A L oimi 
in Eq.Q and choosing the digital recombination gains according to Oii nl allows us to 
obtain a statistics for the quadratures measurement as if a network of beam-splitters 
generating a linear cluster state had been applied to the input squeezed modes. We 
have checked that a numerical solution exists for other choices of the dephasing matrix, 
namely for the one corresponding to the physical dephasings of the cavity modes in the 
experiment of Ref . [15] . 



4. Simple examples of quantum computations 

We now address the emulation of a simple example of quantum computation in the 
measurement-based model (MBQC). Since we are performing homodyne detection on 
all the modes, we are restricted to gaussian operations This is a particularly simple 
framework to start with, because the implementation of gaussian operations in the 
MBQC model does not require an effective adaptivity - usually necessary to implement 
operations deterministically - and the measurements can be performed on all the modes 
simultaneously. The outcomes are recorded, and the "adaptation" of the measurement 
basis depending on the outcomes of previous measurement can be recovered in a further 
classical processing step (e.g., by electronically correcting the photocurrents, see dashed 
box in Figjl]) [B]. 

We consider here single-mode operations. A quantum computation implies first 
coupling an initial quantum state to be processed to a cluster state, and then suitable 
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Figure 4. Fourier transform of an input state through measurement-based quantum 
computation. The input mode is coupled via a beam-splitter (dashed line) to the 
first mode of a three-mode cluster state, and then these two modes are simultaneously 
measured, which teleport the input state onto the second mode. Then, a suitable 
quadrature measurement is performed on the second mode of the cluster, projecting 
the last mode onto the desired output. With the notation used in this picture [TH], the 
angles appearing on each mode specify which quadrature g'jSin^ +p'jCOs£?i shall be 
measured on it in order to implement the desired operation. 

measurements of the quadratures on each mode of the cluster (Fig|4]). The measurement 
of the last mode is not part of the computation but provides a read-out of the output 
state, the "result" of the computation. The coupling between the initial state and 
the first mode of the cluster can be realized by a generalized teleportation [20]: these 

a'i 



two modes are coupled by a beam-splitter interaction, eg. 



in 






with t/gg' 1 = ^ ^ | )• Measuring the quadratures g'^sin^ + p' in cos 9 in and 

g^sin^i + p^cos^i allows to transfer the input state on the second mode of the 
cluster state modulo a transformation M te i e {din,di) controllable via the choice of the 
measurement angles 9 in ,9i |2Uj|| This corresponds to the effective transformation of 
the quadratures 

( q \ ( cos9+ 

= M te le(0in, #l) I „ J With M te i e (9 in , Ox) — — I s i n e°+^Tn6» + 

and 6± = 9 in ±6i (the difference with Eq.(8) of Ref. (20] being due to a different definition 
of the angles 9i n , 0i). In particular, the choice 8 in = 9i = tt/2 teleports exactly the input 
state onto the second mode of the cluster, i.e. 

M«* = f.* = 5) = (i ?)■ (12) 

Let us consider as an example of quantum computation the Fourier transform of the 
initial state. Together with the quadrature displacement D l q (s) = e ls ^ and the shear 
D 2 , q (s) = e isil this is one of the elementary operations for universal one-mode gaussian 
quantum computation [B], and it has been recently implemented experimentally with 
the use of a beam splitter network generating a four-mode cluster state (8]. The Fourier 
transform that we wish to implement should act on the quadratures of the input state 

% The teleported state will be actually display a quadrature shift according to some by-product operator 
which depends on the outcomes of the quadrature measurements on modes in, 1. This can be corrected 
with classical post-processing at the end of the computation. 
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\ J by rotating them, i.e. F = ^ ^ ^ J . For a general cluster state it has been 
shown [6j |20] that measuring the quadrature 

Di2{ s )pDq,2{s) = p + sq = g(xsin9 + pcosO) (13) 

on one of its modes, with g = y/l + s 2 and 9 = arctans, effectively transfers to the 
adjacent (say, right hand) mode a state transformed according to the matrix 



M{s)=\ , I, (14) 




apart from "by-product" displacements which are not accounted for by Eq.(14). 



These may be corrected in the end of the computation, although this is not strictly 



necessary [Sj. Note that the measurement of the quadrature in Eq.(13) corresponds to 
the measurement of the quadrature p after the transformation a — > e t9 a (modulo the 
scale factor g). 

We see that in the simple case under consideration a single computational 
measurement, to be performed on the mode 2 of the cluster state in Fig. [4j is necessary 



after the teleportation step. Indeed, with s 2 = we obtain from Eq.(14) 



M ( S * = °) = K n » ( 15 ) 




leading to M(s 2 = 0)M te i e (6 in = f , X = f) = F. As expressed by Eq.g3f, this 
corresponds to the measurement of the quadrature p on the mode 2. As well as for 
the measurements implementing the teleportation, the corresponding outcome should 
be recorded to implement the final correction step above mentioned [6l [20] . 

Hence, to summarize, in order to implement a Fourier transform on an initial 
state, we have to couple it via a beam-splitter interaction on a three-mode cluster 
state, and then to measure on each mode the quadrature specified by the angles 
(9 i n , 0i, 6*2, 6*3) = (f,f,0, #3) respectively (see Fig. [4]). The measurement on the third 
mode of the cluster yields an outcome probing the implementation of the desired 
transformation, and hence it can be performed according to any quadrature. 

As anticipated, we take the perspective of measuring (simultaneously) the 
quadrature p on all the output modes after suitable rotations have been applied, here 
expressed by the matrix D meas = diag(e l6,in , e 1 ® 1 , e l6>2 , e %dz ) = (i, i, 1, e^ 3 ). Hence the 
unitary transformation which transforms the three p-squeezed modes plus the input 
state (a in ,ai,a 2 ,a 3 ) into the three-mode cluster state coupled with the initial state by 
a beam-splitter interaction reads 

U ti = Ancas ® ^' 3 )(^ 1 ® t/ll n 2 ' 3 ) = Ancas^BS^ 2 ' 3 (16) 

12 3 

where for simplicity we have indicated with the same notation the matrix ' , which 
builds the three-mode linear cluster state from three (infinitely) p— squeezed modes, and 
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i 12 3 

the same matrix times the identity on the fourth mode I <E> ' . 



can be chosen as (see Appendix B ) 
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which, setting 6* 3 = 0, results in the "desired" matrix 
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In order to be able to implement the matrix U t { by means of the multi-pixel homodyne 
detection method, we have to find a combination of A L o and O which satisfies Eq.|4]) 
with U t h = U t { [§} We may assume here appropriate dephasing between the modes, e.g. 



i OPO 



diag(l, 1, 



, i). With the use of Eq.([6]) we obtain 



Among the 2 N = 16 possible solutions for ALOtf an d the corresponding orthogonal 
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matrix O t f satisfying Eqs. wo\ \7y for 
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U t { we can chose for instance the set 
/ 0.89 + 0.46i 
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(20) 
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.(21) 



where we have set ( = arctan a/2. The corresponding orthogonal matrix is 

/ cos (3 — sin j3 cos (3 — sin (3 \ / 

1 sin /3 cos /3 — sin f3 — cos f3 

,J2 — sin (3 — cos j3 — sin j3 — cos f3 

\ — cos (3 sin j3 cos (3 — sin f3 ) 

Note that all the modes are simultaneously measured according to the procedure 
described in the previous section. It is readily verified that O^Alci^G = U%t- 
Hence, with the MPHD we can emulate the formation of a three-mode cluster state, 
the interaction of its first mode with the input mode via a beam-splitter, and the 
measurement of the suitable quadratures to propagate along the cluster the quantum 
computation. Indeed, shaping the phase of the local oscillator on the pixel modes as 
prescribed by the matrix ALOtfi m Eq.([2j and choosing the digital recombination gains 
according to O t n allows us to obtain a statistics for the quadrature measurement of mode 

§ Furthermore, we have to set the global phase of the local oscillator corresponding to the measurement 
of the p quadrature on all the output modes, see Appendix A 
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4 as if this would contain the Fourier transform of the state in the first mode (here in 
a p squeezed state), modulo reinterpreting these result keeping into account the by- 
product operators. Experimentally, addressing the variety of possible input states could 
be realized via proper seeding of the multimode cavity with spatially and temporally 
shaped light. 

The same procedure also allows to implement the quadrature displacement e^ s . 
This is not a symplectic operation, therefore the formalism of Ref. [20J cannot be 
used to determine the angles to measure. In reference to the scheme of FigJIJ after 
teleporting the input state onto the second mode of the cluster state, the measurement 
of the quadrature e~ lh ^'pe lh ^> on this mode allows to implement the gate e~ lh ^\ where 
h(q) is an arbitrary function of the quadrature q [6]. Then, to implement the gate e ts ^ 
the quadrature p S g = e~ ls ^pe lsq = p + s has to be measured. This can be achieved 
by measuring p and adding s to the result. Hence, everything goes exactly as for the 
Fourier transform in terms of the quadratures to measure, except that we have to add 
"s" to the measurement of the quadrature p on mode 2 of the cluster jjj] 

We have also checked that the other elementary operations of the gaussian universal 
set {e^ s ,e^ S ,F,C Z }, namely e^ \ C z may be implemented approximatively, i.e. a 
solution exists for matrices {7th close (in the sense of the matrix distance d{Mi,M<i) = 
1 1 Mi— M2I \f, where | |M| \f = ■sjYlij \Mij\ 2 is the Frobenius norm) to the ones associated 
with these operations. The physical meaning of these approximate operations deserves 
however a detailed study and will be addressed elsewhere. 



5. Conclusion and prospective views 

In this work we have demonstrated the possibility of measuring cluster states and 
implementing gaussian quantum computation in an extremely compact fashion. The 
method is based on the simultaneous measurement of all the (highly) squeezed optical 
modes in a cavity by multi-pixel homodyne detection, and on the classical post- 
processing of the acquired signals, which are multiplied by suitable electronic gains. This 
procedure requires the determination of the suitable phase shape of the local oscillator 
to be employed in the MPHD, as well as the gains to apply to the traces recorded in each 
mode. In particular, as a first example we have provided the explicit solution in terms 
of these experimental parameters which allows to mimic the formation of a four-mode 
linear cluster state. As a simple example of quantum computation, we have considered 
the implementation of a Fourier transform on an input mode. 

It is important to stress that by the Gottesmann-Knill theorem [6], [21] it is known 
that all the results of manipulations with gaussian elements, such as the squeezed 
states and homodyne detection discussed in this work, can be efficiently simulated 

|| Note that the output state obtained by this procedure is |"0 O ut) = Fe l ^ s \ipi n ) (apart from the by- 
product operators arising from the teleportation step). As can be noticed, a Fourier transform multiplies 
the gate e JS * [6J 20j, to be interpreted here as a by-product operator, while in the previous section it 
was considered the desired operation to implement. 
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with a classical computer. In order to perform operations which outperform classical 
processing capability at least a non-gaussian operation is necessary. The inclusion in 
the experimental apparatus of a photon counter which detects the number of photons 
in one mode [22] could lead to the extension of the class of the accessible operations, 
yielding to the implementation of non-gaussian unitaries. 

Finally, we remark that in our set-up the state preparation and the computation 
are performed in a single step, yielding hence a quantum depth equal to 1 - the quantum 
depth being indeed defined as the number of measurement steps that a computation 
task requires, when trying to do as many measurements as possible at the same time. 
This may provide an advantage with respect to the circuit model for performing the 
same operation [23J. 
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Appendix A. Modelization of the operations performed on the modes 

Appendix A.l. Basis and notations 

We introduce here the notations and conventions that we are going to use. The local 
bosonic operator a(p) satisfies 

[a{p),a\p')] = 8{p-p') (A.l) 

and is related to the local quadrature operators by 

q(p) = a\p) + a(p) 

p{p) = i (a\p) - a{p)) . (A.2) 

so that [q(p),p(p')] = 2i5(p — p'). Given an orthonormal basis of wave functions Uk(p) 
satisfying the orthogonality relation 

dpu* k {p)ui{p) = 5 k< i (A. 3) 

the local bosonic operator can be decomposed as 

oo 

a{p) = ^2a Uk u k (p) (A.4) 

fc=0 

with 

a Uk = / dpa(p)u* k (p). (A.5) 



Appendix A.2. Multimode homodyne detection 

We precede now with the modelization of the change of basis which the optical modes 
undergo as an effect of the multi-pixel homodyne detection, as reported in Ref.[I7j. The 
beams splitter mixes the input modes with the local oscillator, yielding the modes 

1 f , . V 

a A,u k - ( a LOJ k + CtuJ 

a B ,u k = ^ (a LOJk - a Uh ) (A.6) 



with 



k 

a A ,u k = / dpu* k (p)a A (p). (A.7) 



Then, light is propagated to the pixels. Let us introduce the output intensities on each 
pixel 

iB,i= / dpa A B (p)a B (p). (A. 8) 



iA,i = / dpa\(p)a A (p) 
J pes, 
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Taking the differences of the intensity signals between couples of corresponding pixels 
gives 

Si — lA,i 1-B,i 

= \Y1 [( a LO,« fc + oi k )(oLO,«, + <) - (°lo, Uk ~ ai k )(a LO , Ul ~ a Ul ) 

k,l 



x / dpu* k (p)ui(p) 

k,l J P> 



dpu* k {p)ui(p) 



(A.9) 



Now we assume that the local oscillator is in an intense coherent state in a certain 
mode that we will indicate by aLO,« LO ( u lo is indeed one of the modes of the basis 
{ui}), and we keep only the terms implying such a mode. Furthermore, we use 



a 



t 

LO,ULO 



(a[ ULO ) = tt 0' being a 3> (d'ljVk. With this we obtain 



a LO,u LO au l 



pes. 



dpu* h0 (p)ui(p) + XX fe aL o,«Lo / d P u l(p) u ho(p) 



P&St 



- y2( a o^u k / dpu* LO (p)u k (p) + a ai k / dpu* k (p)u LO (p) ) . (A.10) 
k V ^pes t ./pes, / 

According to the definition in Eq.(A.5) and to Eq.Q we can write 



dpd(p)v*(p) = Ki I dpd(p)u* LO (p). 

P&Sr 



(A.ll) 



Substituting Eq.( |A.4[ ) in Eq.(A.ll) we obtain 

= / d P u to(p) u j(p) d u 3 = ^2u TiJ d U] , 

where we have defined the matrix realizing the change of bases 

U T i,j = Ki dpu* LO {p)uj{p). 
J pes, 

Remark that the matrix Ut is not necessarily square nor has finite dimension. In terms 
of the pixel modes defined in Eq.Q Equation (A. 10) becomes 

a'*d Vi + a' al 



(A.12) 



(A.13) 



*ol I (e^fo, + iPv t ) + e-^{q vt - ip Vi )) 



= \a' \{q Vi cos(f -p Vi sin^ ) = K|g£°, (A.14) 

where we have set a' Q = \ao\e~ llpo / ' k. We see that the measurement of the quadrature p 
on each output mode requires setting the phase of the local oscillator to <po = 3ti/2. 



Appendix A. 3. Signal recombination 

The pixel modes can be digitally multiplied by real gains and recombined, leading 
to the modes v i (p) = £L g^v'^p) . Since the output functions v i (p) should form an 
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orthonormal basis, we are restricted to gains which are elements of an orthogonal matrix, 
i.e. Oj 



g % y Then the corresponding annihilation operators are 

and we can obtain information on a quadrature as 

s- = 0i >f s 'j = [ a 'oYl 0i iH + a 'o Yl 0i i a \ ) 

3 \ 3 3 J 

a£a v » + a' aj, 



a 



av 



« cos^o -p v " amcpo) 

i i 

u > cos(v9 + <pi) - p v > sin(y? + (p^) 



(A.15) 



(A.16) 



Appendix B. Three-mode linear cluster state matrix 

In this appendix we derive the matrix U^' which builds-up the 3-component linear 
cluster state. For a general cluster state identified by an adjacency matrix V the 
corresponding unitary matrix which brings from the squeezed modes to the given cluster 
state can be chosen as any possible matrix satisfying [21] 

U = X + iY; 
Y - VX = 
XX T + YY T = X; 
X T Y = Y T X; XY T = YX T . 

Setting A = XX T the previous equations leads to the linear system 

VAV = X - A. 

For a three-mode linear cluster state, the adjacency matrix is given by 

/ 1 \ 
1 1 



(B.l) 
(B.2) 



V 



1 



(B.3) 



The linear system in Eq.(B.2) yields in this case the solution 
/ 



A 



2 
3 





V 





1 

3 







2 
3 



(B.4) 



/ 



From Eq.(B.4) we easily obtain the symmetric solution X s such that X* = A, and the 



corresponding U s and V s satisfying Eq.(B.l) 

/ (3+y^) 

X. = 



V 





'-3HV3) 




1 





6 



(3+^3) 

6 



( 





1 

V3 





\ 




1 





1 






V3 


V3 




\ 





1 

v/3 





/ 



; U 8 



( (3+V3) ^_ 

6 ^3 

i 1 

V3 V3 

6 



-3+n/3") \ 
/ 



V3 
(3+^ 
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From this symmetric solution, all the other solutions can be obtained. Indeed, it is 
easy to see from Eq.(B.l) that if X s is a solution, also X = X s O is a solution for any 

12 3 

real orthogonal matrix O. Hence, any matrix U lin ' = X + iY = X s O + iVX s O = 
(I + iV)X s O = U s O can be used to build a cluster state. In order to explicit the free 
parameters induced by this rotation, we will use the Euler parameterization 



O 



' cos ip 
— sin if) 

V 



o 



sin^ 
cos^ 








( 



\ 



cos 9 


— sin 9 



sin# \ 


cos 9 



( 



\ 



cos<p 
— sin( 




sin c 

COS( 








1 



.(B.5) 



With this we obtain U^' 3 



U s O. An example is provided in Eq.(17) of the main text. 



